function [Ftd Ftv] = forceCalculator(pV,pD,pVDot,pDDot,N,kh,kv,lh,lv,b)

%Forces Exterted by the springs
Fv(1:N-1,:)=fSpring(pV(2:N,:),pV(1:N-1,:),kh,lh); %Ventral(N)-Ventral(N+1)
Fd(1:N-1,:)=fSpring(pD(2:N,:),pD(1:N-1,:),kh,lh); %Dorsal(N)-Dorsal(N+1)
Fdv(1:N,:)=fSpring(pV(1:N,:),pD(1:N,:),kv,lv); %Ventral(N).Dorsal(N)	
		
	%Forces acting on Boundary Nodes
Ftd(1,:) = Fd(1,:)+Fdv(1,:); % Dorsal-left BC
Ftd(N,:) = -Fd(N-1,:)+Fdv(N,:); % Dorsal-Right BC
Ftv(1,:) = Fv(1,:)-Fdv(1,:); % Ventral-Left BC
Ftv(N,:) = -Fv(N-1,:)-Fdv(N,:); % Ventral-Left BC
	
i=2:N-1;	
	
% Internal Nodes	
Ftd(i,:) = -Fd(i-1,:)+Fd(i,:)+Fdv(i,:);
Ftv(i,:) = -Fv(i-1,:)+Fv(i,:)-Fdv(i,:);

Ftd = Ftd - b*pDDot;
Ftv = Ftv - b*pVDot;